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ABSTRACT 

We determine the low-redshift field galaxy stellar mass function (GSMF) using an area of 
143 deg^ from the first three years of the Galaxy And Mass Assembly (GAMA) survey. The 
magnitude limits of this redshift survey are r < 19.4 mag over two thirds and 19.8 mag over 
one third of the area. The GSMF is determined from a sample of 5210 galaxies using a density- 
corrected maximum volume method. This efficiently overcomes the issue of fluctuations in the 
number density versus redshift. With Hq = 70kms~^ Mpc~^, the GSMF is well described 
between 10^ and lO^'^M© using a double Schechter function with M* = lO^^ '^'^Mo, 0* = 
3.96 X IQ-^ Mpc"^ ai = -0.35, 0^ = 0.79 x 10"^ Mpc"^ and a2 = -1.47. This result 
is more robust to uncertainties in the flow-model corrected redshifts than from the shallower 
Sloan Digital Sky Survey main sample (r < 17.8 mag). The upturn in the GSMF is also seen 
directly in the z-band and A'-band galaxy luminosity functions. Accurately measuring the 
GSMF below lO^M© is possible within the GAMA survey volume but as expected requires 
deeper imaging data to address the contribution from low surface-brightness galaxies. 

Key words: galaxies: distances and redshifts — galaxies: fundamental parameters — galax- 
ies: luminosity function, mass function 



1 INTRODUCTION 

The galaxy luminosity function (GLF) is a fundamental mea- 
surement that constrains how the Universe's baryonic resources 
are distributed with galaxy mass. Before the advent of CCDs 
and near-IR arrays, the GLF had been primarily measured in 
the B-band jpeltenl [19771: iBinggeli, Sandage. & Tamnianij|l988l : 
lEfstathiou. EUis. & Petersonll 19881 : iLovedav et alj|l992h . More re- 
cently the low-redshift GLF has been measured using thou- 



sands of galaxies in the redder visible bands terown et 
200 ll: iBlanton et all l2003bi) and the near-IR ( Coleetal. 200 IJ 



in 



Kochanek et alj l200lir which more closely follows that of the 



underlying galaxy stellar mass function (GSMF). Furthermore 
the increased availability of multi-wavelength data and spec- 
tra enables stellar masses of gal a xies to be estimated us- 



ing colours ^Larson & TinslevI 1 19781 : Ijablonka & Arimoto 



Bell & de Jond 1200 ih or spectral fitting ([Kauffmarm et al 



Panter. Heavens. & Jimenea 1 20041 : iGallazzi et alj 2005h . and ei 



1992 



2003 
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ther of these method s allow the GSMF to be computed 
dSalucci & Persid Il999l: iBalogh et alj l200lh ISell et alj l2003l : 
iBaldrv. Glazebrook. & Drivej2 008^. hereafter. BGD08). 

Measurement of the GSMF has now become a stan- 
dard tool to gain insights into galaxy evolution with con- 
siderable effort to extend analys es to z ~ 1 i Drorv et al 
l2009l : |Pozzetti etalj 2Q10: Gilbank e taDl201 iLlVulcani et al.l201 1[) 
and higher jEl sneT. Feulner. & Hopp 2008; Kajisawa et al. 2009; 



Marchesini et al. 2009 ; Ca puti et al. ,2011 ; Gonzalez et al., 2011.; 



Mortlock et a D l201lh . Overall the cosmic stellar mass density is 
observed to grow by a factor of 10 betwe en z ^ 2-3 and z = Q 
dPickinson et alJl2003l ; lElsner et al]|2008h with significantly less 
re lative growth in massive > IQ^^M^p, galaxies since z ^ 1- 
2 dWake et alJliooel ; IPozzetti et alj|2010l ; Icaputi et alJl2oTlh . The 
evolution in the GSMF is uncertain, however, with some authors 
suggesting there could be evolution in the stellar initial mass func- 
tion (IMF) (Dave 2008; Wilkins et al. 2008; van Dokkum 2 0o3) 
and considering the rang e of uncertaintie s associated with estimat- 
ing s tellar masses tMaraston et alj|2006l ; IConrov. Gunn. & White! 
I2OO9I) . 

The observed GSMF, defined as the number density of galax- 
ies per logarithmic mass bin, has a declining distribution with 
mass with a sharp cuto ff or break at high masses often fitted 
with a ISchechtej Il97i) function. At low redshift, the character- 
istic mass of the Schechte r break has been determin ed to be be- 
tween 10^° '^ and 10" M0 jPanter et alj2007l ; BGD08; lLi & White! 
!2009 '). The GSMF shape, however, is not well represented by a sin- 
gle Schechter function with a steepening below 10^° Mp giving rise 
to a dou ble Schechter function shape overall (BGD08). !Peng et"an 
(l2010bl) note that this shape arises naturally in a model with simple 
empirical laws for quenching of star formation in galaxies. This is 
one example of the potential for insights that can be obtained by 
studying the inferred GSMF as opposed to comparing observations 
and theoretical predictions of the GLF; though we note that it is 
in some sense more natural for theory to predict the GLF because 
model galaxies have a 'known' star-formation history. 

Abundance matching between a theoretical galactic halo 
mass function and a GLF or GSMF demonstrates that, in or- 
der to explain the GSMF shape, the fraction of baryonic mass 
converted to stars increases with mass to a peak before de- 
creasing ( Marinoni & Hudson 2002; Shankar et al. 2006; B GD08 ; 
!Conrov & Wechsler. .200a ; iGuo et ah ,20 IQ ; .Moster et all hOld) . 
Galaxy formation theory must explain the preferred mass for 
star formation efficiency, the shallow low-mass end slope com- 
pared to the halo mass function, and the exponential cutoff at 
high masses dOppenheimer et ^ !2010l) . At high masses, feedback 
from active galactic nuclei has been invoked to prevent cooling of 
gas leading to sta r formation (B estetal. 200i !Keres et aD !2005l ; 
!BoweretalJ!2006! ; !Croton et alj|2006() . The preferred mass scale 
may correspond to a halo mass of ^ 10^^ above which gas 
becomes more readily shock heated jPekel & Bimboiml2006l) . To- 
ward lower masses, supemovae feedback creating galactic winds is 
thoug ht to play a rnajor role in regulating star formation jLarsonI 
'1- Dekel & Silk 'l98^; !Lacev & Sillj !l99ll; !Kav et al.! BoO^ 



11974 ; 

while iOppenheimer et alj have argued that it is re-accretion of these 
winds that is critical in shaping the GSMF. Others have argued that 
star formation in the lowest mass haloes is also suppressed by p ho- 
toionization (Efstathiou 1992; Thoul & Weinberg 1996; Somervill^ 
l2002i) . in particular, to expl ain the number of satellites in the Local 
Group teenson et al]!2002h. 

Recently, !Guo et a\. I1 I20T l|) used a s emi-analytical mode l ap- 
plied to the Millenium Simulation (MS) jSpringel et alJ!2005!) and 



the higher resolution MS-II teovlan-Kolchin et alj[2009h to pre- 
dict the cosmic-average 'field' GSMF down to lO^M©. They find 
that the GSMF continues to rise to low masses reaching > 0.1 
galaxies per Mpc'' for lO'^-lO^M©, however, they caution that their 
model produces a larger passive fraction than is observed amongst 
the low-mass population. iMamon et al. il) apply a simple one- 
equation prescription, on top of a halo merger tree, that requires star 
formation to occur within a minimum mass set by the temperature 
of the inter-galactic medium. Their results give a rising baryonic 
mass function down to lO^'^M©, with the peak of the mass func- 
tion of the star-forming galaxy population at ~ lO^M©. We note 
that it is useful for theorists to predict the field GSMF of the star- 
forming population because measuring this to low masses, while 
challenging, is significantly easier than for the passive population. 

Measurements of the GLF reaching l ow luminosi- 
ties have been made for the Local Group ('Koposov et alj 
2008h. 



selected groups ([jr entham & Tully 2002 



Chiboucas. Karachentsev. & Tullvl |2009|) . clusters dSabatini et alJ 
2003 ; !Rines & Gelled !2008h and superclusters jMercurio eTo] 



20061) . To accurately measure the cosmic-average GSMF, it is 
necessary to survey random volumes primarily beyond ~ lOMpc 
because: (i) at smaller distances, the measurement is limited 
in accuracy by systematic uncertaint ies in distances to galaxies 
( iMasters. Havnes. & Giovanelli! !2004l) ; and (ii) a local volume 
survey is significantly biased, e.g., the B-band luminosity density 
out to 5 Mpc is about a factor of 5 times the cosmic average 
(using data f r om Karachentsev et al] 20041 in comparison with 
!Norberg et all !2002! ; !Blanton et alj !2003bl) . The Sloan Digital 
Sky Survey (SDSS) made a significant breakthrough with red- 
shifts obtained to r < 17.8 mag and multi-colour photometry 
("Stoughton et al.' 2002"), in particular, with the low-redshift sample 
described by B Ian ton et al . (2005a3. In order to extend and check 
the low-mass GSMF (BGD08), it is necessary to go deeper over a 
still significant area of the sky. 

Here we report on the preliminary analysis to determine the 
z < 0.06 GSMF from the Galaxy And Mass Assembly (GAMA) 
survey, which has obtained redshifts to r < 19.8 mag currently 
targeted using SDSS imaging but which ultimately will be updated 
with deeper imaging. The plan of the paper is as follows. In § |2] 
the data, sample selection and methods are described; in §|3] the 
GLF and GSMF results are presented and discussed. Summary and 
conclusions are given in §|4| 

Magnitude s are corrected for Galactic extincti on using the 
dusts maps of ISchlegel. Finkbeiner. & Davis! ( !l998h , and are k- 
corrected to compute res t-frame colours and absolute magnitudes 
using KCORRECT v4_2 telanton et~an!2003a! ; !Blanton & Roweij 
! 20071) . We assume a flat ACDM cos mology w i th H p = 
70 km Mp c~^ and fl^.p = 0.3. The !Chabriei] ( !2003!) IMF 
(similar to the i Kroupa! !200l! IMF) is assumed for stellar mass 
estimates. Solar absolute magnitudes are taken from table 1 of 
!HilletalJ i l2oToh . and mass-to-light ratios are given in solar units. 



2 DATA AND METHODS 

2.1 Galaxy And Mass Assembly survey 

The GAMA survey aims to provide redshifts and multi-wavelength 
images of > 250 000 galaxies over > 250 des^ to r < 19.8 mag 
toriver et al.l2009l!201l! ; !Baldrv et alj2oTol) . A core component of 
this programme is a ga laxy redshift surve y using the upgraded 2dF 
instrument AAOmega jSharp et alj|2006i) on the Anglo- Australian 
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Telescope. The first three yea rs of the redshift survey have been 
completed toriveretal]|201lh and these data are used here. The 
target selection was over three 48 deg^ fields centred at 9 h (G09), 
12 h (G12) and 14.5 h (G15) on the celestial equator. The limiting 
magnitudes of the main survey were r < 19. 4 in G09 and G15 , 
r < 19.8 in G12, z < 18.2 and Kab < 17.6 dSaldrv et alj|2010h . 
It is only the r-band selection that is used in the current analysis 
because the near-IR selections add mainly higher-redshift galaxies. 
Each area in the survey was effectivel y 'tiled' 5-10 times wi th a 
strategy aiming for h igh completene ss ( Rqbotham et ai] |2010l) . In 
other GAMA paper s, Lovedav et alj ( l2012 h determines the ugriz 
GLFs, [Driver etaL I i2012l) determines the cosmic s pectral energy 
distribution from the far ultraviolet to infrared, while lSrough et aU 
looks at the properties of galaxies at the faint end of the Ha 

GLF. 

The target selection was based primarily on SDSS DR6 
jAdelman-McCarthv et alj |2008[) with the A'-ban d selection us- 
ing U KIRT Infrared Deep Sky Survey (UKIDSS, Lawrence etal] 
20071). and star- galaxy separation using both surveys (see 



Baldrv et alJlioiol for details). Quality control of the imaging se- 
lection was done prior to the redshift survey to remove obvious 
artifacts and deblended parts of galaxies. An update to the target 
list was performe d to remove target s with erroneous selection pho- 
tometry (§ 2.9 of lDriveretZll201ll) . For this paper, further visual 
inspection was made of low-redshift 'pairs' with measured veloc- 
ity differences < 300 kms~^. This resulted in ~ 50 objects being 
reclassified as a deblended part of a galaxy. Further inspection was 
made of targets with faint fibre magnitudes, reclassifying ~ 100 
objects as not a target. After this, the r-band magnitude limited 
main survey consists of 1 14 360 targets Q 

Various photometric measurements are used in this paper. 
The selection m agnitudes are SDSS Pe trosian magnitudes from the 
PHOTO pipeline jStoughton et alj2002l) . In order to obtain matched 
aperture photometry from u- to K-h&nA, the imaging data from 
SPS S and UKIDSS were repr ocessed and run through S EXTRAC- 
TOR teertin& Amoutill99d) . The details aie given in Hill et alj 

l|) and here we use the r-defined AtJTO magnitudes primarily 
for colours: these use elliptical apertures. SEXTRACTOR fails to 
locate some genuine sources that were identified by PHOTO, how- 
ever, and for these we use Petrosian colours. Finally, an estimate of 
total luminosity is obtained using Sersic fits extrapolated to IQRe 
(ten times the half-light radius). This p rocedure uses a fe w soft- 
ware packages inclu ding GALFIT v er. 3 jPeng et al.ll2oToah and is 
described in detail bv lKelvin etall feOllh . 

Spectra for the GAMA survey are taken with the AAOmega 
spectrograph on the Anglo-Australian Telescope (AAT), coupled 
with various other public survey data and some redshifts from the 
Liverpool Telescope. The AAT d ata were reduced using 2dFdr 
jCro om. Saunders. & Heald 2004) and the redshifts determined us- 
ing RUNZ (Saunders, Cannon, & Sutherland 2 004) . The recovered 
redshift for each spectrum is assigned a quality Q from 1 (no red- 
shift) to 4 (reliable). These are later updated based on a compara- 
tive analysis between different opinions of a large subset of spectra. 



^ The sample was derived from the GAMA database table 
TilingCatvie with SURVEY.CLASS J: 6. GAMA AUTO and 
Sersic photometry were taken from tables rDef PhotomvOl and 
SersicCatv07; stellar masses from StellarMassesv03; red- 
shifts, qualities and probabilities from SpecAllvOS; flow-corrected 
redshifts from DistancesFramesvO 6; and photometric redshifts 
from Photozv3. SDSS photometry was taken from SDSS table 
dr6 .PhotoObj. 



From this process, a best redshift estimate and the probability {pz) 
of whether this is correct are assigned to each spectrum. The new Q 
values are based on these probabilities (formally called nQ, § 2.5 
of Driver et al. 201 1). Where there is more than one spectrum for 
a source, the redshift is taken from the spectrum with the highest 
Q value. For the r-band limited main sample, 93.1%, 3.0% and 
3.4% have Q ^ 4, Q = 3 and Q = 2, respectively. In general, 
redshifts with Q ^ 3 are used, however, Q = 2 can be considered 
when there is agreement with a second spectrum that was measured 
independently or when there is reasonable agreement with an inde- 
pendent photometric redshift estimate. 



2.2 Stellar mass estimates 

Stellar masses were computed for GAMA t argets using the ob- 
served AUTO matched aperture photometry jffill etalj|20Tll) for 
the ugriz bands. These were fitted using a grid of synthetic spec- 
tra wi th exponentially declining star formation hi stories produced 
using lBruzual & Chariot (2003) models with the IChabrieJ b003h 
IMIrl and the lCalzetti et al. (2000) dust obscuration law. The stellar 
masses were determined from probability weighted i ntegrals over 
the full range of possibilities provided by the grid. See lTavlor et alj 
(201 1) for details of the method. For the fitted stellar masses in this 
paper, we use the stellar mass-to-light ratios (M/L) in i-band ap- 
plied to the j-band Sersic fluxes. Where M/L values are not avail- 
able (2 % of the low - redsh ift sample), we use the colour-based rela- 
tion of lTavlor et al. U20Ilh to estimate M/Li . 

The 95% range in M/L,; from the fitting bv lTavlor et al] i201lh 
is 0.5-2.0 (M0/L0) for high luminosity galaxies (Li > 1O^°L0) 
and 0.2-1.6 for lower luminosity galaxies (z < 0.06). The net un- 
certainty on an individual stellar mass e stimate can be large, e.g., a 
factor of two or 0.3 dex as estimated bv lConrov e t al. ( 2009). Note 
though that the impact of uncertainties is more important when con- 
sidering evolution in the GSMF than when considering the shape 
of the GSMF at a single epoch as in this paper. The latter primarily 
depend s on the diiferential systematic uncertainty between popu- 
lations. I^^loretal] |201o2) estimated that the net differential bias 
was < 0.12 dex based on comparing stellar and dynamical mass 
estimates. The change in M/Li between red and blue galaxy pop- 
ulations can be approximated by a colour-based M/L relation. The 
effect of changing the slope of this relation is considered in §[331 

We note that the reason that M/L correlates with colour at all 
well is that the M/L of a stellar population increases as a popula- 
tion reddens with age or dust attenuation. ISell & de Jond j200lh 
noted that errors in dust estimates do not strongly affect stellar 
mass estimates. In other words, the vectors in M/L versus dust red- 



dening run nearly parallel to those determin ed for age reddening . 
[Driver et d]j2007i) . using the dust models of IPopescu et alj j2000l) 



Tuffs et al. H2004h . confirmed this with the largest deviation for 



and 



edge-on systems. 



2.3 Distances 



The GAMA survey, as with most large redshift surveys, provides 
the heliocentric redshift as standard. In many cases it is sufficient 
to assume that this is close to the cosmological redshift. For the 
GAMA regions at low redshift, however, it is not. 



^ Stellar masses derived using the Chabrier IMF are about 0.6 times the 
masses derived assuming the Salpeter IMF from O.IMq to IOOMq. 
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Figure 1. The relation between distance and CMB-frame redsliift from the lTonrv et alj )2000l) flow model. The three GAMA regions are shown in different 
panels. The black line represents the central sight line while the grey region shows the variation over each GAMA region. 



First we convert heliocentric redshifts to the cosmic mi- 
crowave background (CMB) frame: 



1 ^cmb (-^ ^hclio)(l ^~ ^sun,con 



(1) 



where czgun.comp is the component of the Sun's velocity toward 
the object in the CMB frame. We use z^un = 0.001231 (369 km/s) 
for the CMB dipole in the direction of RA = 168.0° and DEC = 
—7.2° jLineweaver et al. 1996). For the GAMA survey, this leads 
to average corrections for the heliocentric velocity (c Zsun.comp) of 
+303, +357 and +236kms"\ in G09, G12 and G15, respec- 
tively. 

In the absence of flow information, the CMB frame red- 
shift is a preferred estimate of the cosmological redshift at z > 
0.03; the velocity of the Local Group (LG) with respect to 
the CMB has bee n attributed to superc lusters at lower redshifts 
jTonrv et alj |2000| : lErdogdu et al] |2006|) . However, it should be 
noted that large-scale bulk flows have bee n claimed by, for ex- 
ample, l^^tyns^^eldmm^^^&Hudso^ ([ioO^. To a ccount for fl ows 
in the nearby Universe, we use the flow model of IXonrv et al.l lin- 
early tapering to the CMB frame between z = 0.02 and z = 0.03. 
Figure [T] shows the relation between the flow-corrected and CMB 
frame velocity at z < 0.01. The main feature is the difference be- 
tween velocities in front of and behind the Virgo Cluster for sight 
lines in G12. For each object sky position, the flow-corrected red- 
shift is obtained by computing the model CMB frame redshifts 
(2cmb,m) for a vcctor of cosmological redshifts {zm from —1000 
to +1000 kms~^ in steps of 10 around the observed Zcmb)- The 
flow-corrected redshift (z) is then given by the weighted mean of 
Zra values with weights 



■ exp 



2cmb ) 



2cr2 



(2) 



where a is taken as 50kms^^. This small value is chosen so that 
the result is nearly equivalent to using the one-to-one solution for 
flow-corrected redshift from 2:cmb, which is mostly available, and it 
cor responds to a typic al redshift uncertainty from the GAMA spec- 
tra jPriver et al.ll201 ih . It is only around 1500 km/s in G12 where 
the method is necessary to provide a smoothly varying weighted 
average between the three solutions. 

Figure [2{ a) shows the difference in distance modulus (DM) 
using the flow-corrected z compared to using Zhoiio versus redshift. 
Note that the correction to a DM can be larger than 0.5 mag; the 
direction of G12 in particular is within 20° of both the CMB dipole 



and the Virgo Cluster (cf. fig. 5 of IJones et al.]|2006l for the South- 
ern sky). The DM uncertainty for each galaxy was estimated by ap- 
plying changes in heliocentric velocity of ±180 km/s and recom- 
puting the flow-corrected distances. T his corresponds to the cos- 
mic thermal velocity dispersion in the iTonrv et al.l j2000l) model, 
i.e., the velocity deviations after accounting for the attractors. Fig- 
ure [2{b) shows the DM uncertainties, which are taken as half the 
difference between the positive and negative changes. The uncer- 
tainty is less than 0.2 mag at 2 > 0.007, while at the lower redshifts 
the uncertainty can be quite large especially in G12 because of the 
triple-valued solution caused by Virgo infall [Fig.fHb)]. 

To show the significance of using different distances, the r- 
band GLF was computed using these flow-corr ected redshifts , he- 
liocentric fra me, LG frame dCourteau & van den Berghlll999l) and 
CMB frame jLineweaver et al.l 1 9961) . Figure [3] shows the 4 resul- 
tant GLFs. The computed number densities are significantly differ- 
ent at magnitudes fainter than ~ —15 because of variations in the 
estimated absolute magnitudes, the sample, and the volumes Vinax- 
The estimates of number density are lower when using the flow- 
corrected redshifts or CMB frame, with respect to heliocentric or 
LG frame, because of larger dis tances and thu s higher luminosities 
with lower 1 /Kna x weighting. Masters et alj ( [200 4) noted that the 
iTonrv et alj 1 I2OOOI) model works well toward Virgo although pos- 
sibly at the expense of the anti- Virgo direction, but in any case, 
the model i s suitable for the GAMA fields. Hereafter, we use the 
iTonrv et al] flow-corrected redshifts. 

2.4 Sample selection 

In addition to the survey selection described in lBaldrvetalJ ( l2010l) . 
the sample selection is as follows: 

(i) rpct < 19.4maginG09orG15,orrpot < 19.8 mag in G12; 

(ii) redshift quality Q ^ 3, Q = 2 when there was agreement 
with a second independent spectrum of the same target (within a 
velocity difference of 450kms~^), or Q = 2 with > 0.7 
and agreement with a photometric redshift estimate [within 0.05 
in 5z/{l + z)]\ 

(iii) 0.002 < z < 0.06 (flow corrected), comoving distances 
from 8.6 Mpc to 253 Mpc; 

(iv) physical Petrosian half-light radius > 100 pc. 

The magnitude limits define the r-band limited main sample 
(114 360) with 98.3% (112393) satisfying the redshift quality cri- 
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z (flow cor'i'ected r'edahift) 

Figure 2. (a): The difference in the distance moduli derived from z and 
^helio- The lines connect the locations of the galaxies for each region, which 
have been sorted by redshift. The vertical dashed line shows the low-redshift 
limit for the further analyses in this paper, (b): An estimate of the uncer- 
tainty in the the distance moduli. The distance moduli were recomputed 
after adjusting the heliocentric redshift by ±180 kms~^ and half this dif- 
ference is plotted. 
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Figure 3. Comparison between galaxy luminosity functions computed us- 
ing different approximations for the cosmological redshift. The samples 
were selected over the redshift range 0.003 to 0.06, and a standard Vmax 
method was used. The solid line with error bars shows the GLF using the 
flow-corrected distances (offset xO.2, also shown with no offset without 
eiTor bars). 



teria. The redshift range reduces the sample to 5217 (50 of these 
were included because of the Q — 2 agreement tests). A further 7 
sources are rejected by the half-light radius criterium giving a pri- 
m ary sample of 5210 g alaxies. The 100 pc lower limit corresponds 
to lGilmore et alj ( l2007h 's division between star clusters and galax- 
ies. However, on inspection the rejected sources were simply as- 
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Figure 4. Absolute magnitude versus redshift. The black and blue points 
represent galaxies in G09 and G15 to r < 19.4 while the red points repre- 
sent galaxies in G12 to r < 19.8. 



signed an incorrect redshift and should be either stellar or at higher 
redshift than our sample limit. 

Figure |4] shows the distribution of the primary sample in AU 
versus redshift. Note not all the redshifts come from the GAMA 
AAOmega campaign with the breakdown as follows: 2671 GAMA, 
2007 SDSS, 444 2dF Galaxy Redshift Survey, 64 Millennium 
Galaxy Catalogue, 10 6dF Galaxy Survey, 6 Updated Zwicky Cat- 
alogue, 6 Liverpool Telescope, and 2 others via the NASA/IPAC 
Extragalactic Database. 

The Petrosian photometry, used for selection of the sample, 
is highly reliable having undergone various visual checks. The ex- 
ception is for overdeblended sources. For these the deblended parts 
have been identified and associated with a target galaxy. The r- 
band Petrosian photometry of these overdeblended sources is re- 
computed by summing the flux from identified parts. About 100 
galaxies have their Petrosian magnitude brightened by > 0.1 mag 
from this, with 14 brightened by more than a magnitude (the target 
part has not been assigned the majority of flux in a few cases). It 
is important to do this prior to calculating Vmax because a nearby 
galaxy that is deblended into parts would not be deblended nearly 
so significantly if placed at higher redshift. 

2.5 Density-corrected maximum volume method 

A standard method to co mpute binned G LFs is through weighting 
each galaxy by 1/Vmax dSchmidtll 19681) , which is the comoving 
volume over which the galaxy could be observed within the survey 
limits (zmax is the corresponding maximum redshift). In the pres- 
ence of large-scale structure, large variations in the number den- 
sity versus redsh ift, this method can distort the shape of the GLF 
tefstathiou et alj 1988). Figure |5] shows the large-scale structure in 
and around the GAMA regions. There are a few substantial over- 
densities and underdensities as a function of redshift within each 
region. 

In order to compute binned GLFs undistorted by radial vari- 
ations in large-scale structure, a density-corrected Vj^ax method is 
used. This is given by 

■^'-^^ Ekl: 



; Ci 



(3) 



where is the completeness factor assigned to a galaxy; and the 
corrected volume is given by 
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Figure 5. Redshift distribution in and around tlie GAMA regions. Only galaxies witli r < 17.8 are sliown corresponding to the SDSS main galaxy sample 
limit (but including redshifts from all surveys). Note the low number of galaxies at low RA, to the left of G09 in this figure, is because of the lack of redshift 
survey coverage from SDSS, 2dFGRS and GAMA. 



pddp[Zl\ Z2) 

where pddp{za',Zb) is the number density of a density-defining 
population (DDP) between redshifts Za and zt, z\ = 0.002 is 
the low-redshift limit, and Z2 = 0.06 is the high-redshi ft limit of 
the sa mple. This method is also described in § 2.7 of Ba ldry et"an 
( I2OO6 ) where the density-corrected volume is given by /vVlnax- 
To calculate this, we first treat G09-I-G15 and G12 separately be- 
cause of the different magnitude limits, except that we use a single 
value for pddp(zi; 22), which is taken to be the average density of 
the DDP over all three regions. The DDP must be a volume-limited 
sample and we use Mr < —17.9 (Fig.|4l(. Figure [6] shows the rel- 
ative number density [fv{z) = pddp{zi; z)/pddp{zi; Z2)] for the 
separate samples. 

The redshift upper limit of 0.06 allows sufficient statistics 
to be obtained at the bright end to fit the knee of the GLF or 
GSMF while at the same time allowing the use of a single DDP 
that can be used to reliably measure Knax for galaxies as faint as 
Mr ~ —13. Raising the redshift limit to 0.1 would improve the 
bright-end statistics at the expense of using a DDP with a limit that 
is 1.1 mag brighter, which is less accurate for determining V^ax 
values. We note, however, that it is possible to use a series of over- 
lapping volume-limited samples to improve the accuracy of Vj^ax 
(e.g. Mahtessian 201 1 'sewed' three samples together). For the pur- 
poses of keeping a simple transparent assumption and mitigating 
against even modest evolution, for this paper we use a single DDP 
with z < 0.06. 

The step-wise maximum likelihood (SWML) 
tefstathiou et alj Il988l) method can also be used to compute 
a binned GLF that is not distorted by large-scale structure 
variations. In fact, computing the density binned radially and 
the binned GLF using a maximum likelihood method can be 
shown to be equivalent to a density-corrected l^niax method (§ 8 
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Figure 6. Variation in number density of the DDP used in the density- 
corrected Vjjjax method. The lines represent the number density of the 
Mr < —17.9 volume-limited sample divided by 0.0148 Mpc^"^ from 
z = 0.002 up to the redshift shown on the X'-axis. 



of lSaundersetai]|l99(]| : IColell201lh . This is reassuring but not 
surprising given that both SWML and density-corrected V^ax 
methods assume that the shape of the GLF remains the same 
between different regions. This is not exactly true but the resulting 
GLF is a weighted radial average. This is seen more transparently 
in the density-corrected V^ax method. The real advantage here 
is that V^niax need only be calculated for each galaxy using the 
selection r-band Petrosian magnitudes after which the GLF 
(or GSMF) can be determined straightforwardly using different 
photometry. When calculating the GLF in a different band (or the 
GSMF) there is no colour bias in a bin unless a population with a 
certain colour is only visible over a reduced range of luminosity 
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Figure 7. Compaiison between standard maximum volumes and corrected 
volumes. The black and grey lines represent Vmax; note the small scatter is 
due to dilferential fc-corrections between z and Zmax- The blue and red lines 
represent V^^ax- ^ minimum value of 100 Mpc^ was set for V^^ax because 
this volume could only be expected to contain typically 1 or 2 galaxies of 
the DDR 



(mass) within the bin. Note also the GAMA DDP sample is highly 
complete, which means that the calculation of pddp is robust. 

Figure |7] shows a comparison between V^^ax ™d Knax- For 
example, note the flattening of Knax in G12 brighter than —15.2 
(red line). This corresponds to the overdensity at 2; ~ 0.022 with 
the underdensity beyond. Brighter galaxies can be seen further but 
the corrected volume rises slower than the standard Knax because 
the DDP is underdense beyond. 

In order to estimate GLFs, the completeness is assumed to 
be unity (d = 1) in this paper with the area of the survey being 
143 deg^ (one third of this for each region). Figure [8] shows the 
i-band GLF computed using the different volume correction meth- 
ods. The Kniax method produces much better agreement between 
the regions than the standard VInax method. The remaining differ- 
ence between the regions, below < lO^L© in particular, may be the 
result of the GLF varying between environments or uncertainties in 
the distances. The grey lines in Fig. [8] represent the GLF using a 
combined volume over all regions. This is obtained by modifying 
Pddp(2i; Zmax.i) Vmax.i in Eq.|4]to be a sum over afl three regions 
for each galaxy with Zmax.i being different in G09-I-G15 (r < 19.4) 
compared to G12 (r < 19.8) (see also Avni & Bahcall 1980 for 
combining samples with different effective volumes). Hereafter, 
this combined V^ax is used. Note also we show GLFs using so- 
lar luminosities because we are working towards the GSMF. 



3 RESULTS AND DISCUSSION 

3.1 Galaxy luminosity functions in the i-band 

The S/N in the i-band is significantly higher than the SDSS z-band 
or any of the UKIDSS bands for galaxies in our sample. Thus we 
use the i-band as the fiducial band from which to apply stellar mass- 
to-light ratios. First we start by looking at the i-band and com- 
paring the GLF taken with different photometric apertures. Fig- 
ure Oa) shows the i-band GLF using photo metry from (i) S DSS 
pipeline PHOTO, (ii) SE XTRACTORa s run b y lHTll et al J jioTH) and 
(iii) GALFI T as run by iKelvin et alj ( [201 l|). For comparison, the 
result from lLovedav et al] ( |2012|) at z < 0.1 using Petrosian mag- 
nitudes and SWML method is also shown (here computed slightly 
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Figure 8. Comparison between the standard Vmax (upper data) and density- 
corrected y^ax (lowsr data) for determining the i-band GLF. The triangles 
show the GLFs from the standard method (offset X 100). Note the major 
discrepancies between the regions. The circles show the GLFs from the 
Vjjjax method. The grey lines (no offset and offset xlOO) represent the 
GLF using a combined V^ax °^er all three regions. 



differently to their paper, with fc-corrections to z = and with no 
completeness corrections). 

The difference between the GLFs in Fig. |9ta) are generally 
small except for the SWML z < 0.1 GLF which is lower around 
the 'knee'. The 2: < 0.1 GAMA volume is known to be underdense 
by about 1 5% with respect to a larger SDSS volume (see fig. 20 of 
i Drive r et"al. 201 1) whereas the z < 0.06 GAMA density is similar 
to the SDSS volume. When the SWML method is applied to a 2 < 
0.06 GAMA sample, there is significantly better agreement with 
the density-corrected Kmax LF as expected. Thus the 2; < 0.1 LF 
has a different normalisation and shape primarily because the i- 
band LF is not exactly universal between different environments. 

The faint end differences in Fig.|9{a) are generally not signifi- 
cant (cf. error bars in Fig.[8]l. At the bright end, the differences are 
because the AUTO apertures and Sersic fits are recovering more flux 
from early-type galaxies than the Petrosian aperture. 

Figure |9jb) compares the GAMA result using PHOTO Pet- 
rosian magnitudes with results using th e SDSS NYU-VAGC low- 
redshift sample (0.0033 < z < 0.05; iBlanton et allbOOSbl) . Ig- 
noring the differences below 1 0^'^Lq, which are be cause of the 
differing magnitude limits, the iBlanton et al] ( l2005al) GLF (DR2) 
gives a higher number density below lO^L©. This can be at least 
partly explained by the distances used . The NYU-VAGC uses dis- 
tances from the IWillick et model, tapering to the LG 
frame beyond 90 Mpc. The black dotted line in Fig. [^b) repre- 
sents the i-band GLF calculated using the method and sample of 
BGD08 (DR4) with the NYU-VAGC distances, while for the black 
dashed line the di stances were changed to those derived from the 
iTonrv et alj ilOOd) model. The latter model gives on average 10%, 
and up to 30%, larger distances at 2 < 0.01. The DR4 result with 
the adjusted distances is in better agreement with the GAMA re- 
sult. Note that GAMA galaxies with luminosities ~ IO^Lq have 
a median redshift of 0.02 compared to 0.006 for the NYU-VAGC 
sample. Thus the GAMA result is less sensitive to the flow model 
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Figure 9. Comparison between i-band GLFs. (a): The red lines, black line and blue dashed line repr esent GLFs computed using the same weights (l/V^a^) 
but different aperture photometry from SDSS pipeline PHOTO, SEXTRACTOR ( thelHill et alj|201ll G AMA photometry) and GALFIT (the lKelvin et al.ll20li] 
Sersic fits). The blue solid line shows the result using the SWML method as perjLoveda^etalJj2012 ) applied to a 2: < 0.1 sample, (b): The red dashed line 
is unchanged. The blue line is from the uncorrected (for completeness) GLF of iBlanton et al.l {20O5i I, and t he black dotted line is computed using the same 
method and sample from BGD08. The black dashed line is the same sample with the distances adjusted to the lTonrv et al.l j2000l) flow model. 
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Figure 10. GAMA data in comparison with a mo del LP. The model L P was 
derived by H. Kim & C. Baugh using GALFORM iBower et al.l2006l) . 



at these luminosities. See lLovedav et al.| ( 1201 2h for more details on 
the GAMA GLFs. 

Figure [TO] shows the GAMA i-band GLF with error bars in 
comparison with the GLF from a semi-analytical model. The latter 
was derived by H. Kim & C. Baugh (private communication) us- 
ing an implementation of GALFORM similar to Bower et al. ( 2006). 
The mass resolution of the halo merger trees was improved and the 
photo-ionisation prescription was changed so that cooling in haloes 
with a circular velocity below 30 km s~ ^ (previously 50 km s~ ^ ) is 
prevented after reionisation (2 = 6). There is reasonable agreement 
between the model and data, however, the model LF is higher par- 
ticular below lO^'^L©. At low luminosities, it is expected that the 
GAMA data are incomplete because of surface brightness issues 
and the LF data points are shown as lower limits (justified in the 
following section, ii l3.2b . 



3.2 Surface brightness limit 

In addition to the explicit magnitude limit, there is an implicit and 
imprecisely-defined surface brightness (SB) limit that plagues mea- 
surements of the fai n t end of a GL F ( Phillipps & Disney 198^; 
ICross&Driverll2002h .l Blanton et al. I (2005 a) estimated the impact 
on the SDSS GLF, determining a completeness of about 0.7 at 
^r,50 = 23.0 mag arcsec"^ where this is the SB within the Pet- 
rosian half-light radius. Three sources of incompleteness were con- 
sidered: photometric incompleteness determined from simulations 
that put fake galaxies in frames run through PHOTO, tiling incom- 
pleteness because some of the SDSS area was targeted on versions 
of PHOTO where the deblender was not performing optimally, and 
spectroscopic incompleteness. The tiling incompleteness is not an 
issue here, the issues associated with the photometric incomplete- 
ness may be less severe at the GAMA faint limit because for a 
given SB the galaxies are smaller, meaning fewer problems with 
deblender shredding and sky subtraction, and the spectroscopic 
incompleteness can be mitigated by repeated observations of the 
same target where necessary. 

Figure [TT] shows the SB versus stellar mas s distribution (with 
masses from the colour-based M/L relation of Taylo r et alj|201ll . 
see ii l3.3t . It is difficult to determine when the input catalogue be- 
comes incomplete. Judging from the slightly higher mean SB at 
lO'^-lO'-'Mo in the GAMA sample compared to SDSS, we expect 
that incompleteness becomes si gnific a nt for surface brightnesses 
slightly fain t er tha n the iBlanton et alj j2005ah estimate. Recently 
iGeller et al] j201ll) analysed a 0.02 < z < 0.1 sample from the 
Smithsonian Hectospec Lensing Survey (R < 20.6). They deter- 
mined a linear relation between SB and Mr for the blue popula- 
tion. This is shown in Fig.[TT]after converting to Ho — 70, assum- 
ing /ir,50 = SB_R,5o + 0.3 and M/Lr — 0.5, which is an aver- 
age value for a star-forming low-mass galaxy. The average GAMA 
SB-mass relation falls below this relation at M < 1O*M0, which 
is where we expect incompleteness to become significant. Rather 
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Figure 11. Surface brightness versus stellar mass. The black diamonds and 
bars represent the GAMA sample used in this paper, while the grey squares 
and bars show the results from an SDSS z < 0.05 sample (as per fig. 4 
of BGD08). The vertical bars represent the measured scatter, ztlcr, around 
the median in bins of 0.5 dex. /ir,50 is the SB within the Petrosian half- 
light radius. The grey dashed-line region represents the expected area of low 
completeness. The dash-and-dotted line shows the blue population relation 
converted from lOeller et alj j201 ih . 

than attempting to correct for this incompleteness, we assume that 
the GSMF values below 10* M© are lower hmits. For the i-band 
LP, values below lO^'^L© were taken to be lower limits (Fig. llOt 
because the M/Li of dwarf galaxies around lO ^Mp) is typically less 
than 0.8 from the fitting of IXavlor et all bOllh . 

3.3 Galaxy stellar mass functions 

Various authors have suggested that M/L in the i-band or 
z-ban d correlates most usef u lly w it h a — i |Gallazzi & Belli 
I2OO9I : IZibetti. Chariot. & Rid l2009l : iTavloret alJ l2010bl) . The 
parametrization is usually linear as follows 

log{M/L,) = a + h{g^i) (5) 

where M is the stellar mass and Li is the luminosity in s olar units. 
Howev er, estimates of a and 6 can vary considerably. iBell et all 
( l2003h gi ve a = -0.152 and h = 0.518 while reading off from 
figure 4 of lZibetti et id] 12009^ gives a ~ — 1 and 6 ~ 1, though the 
latter is for resolved par ts of galaxies. From fitting to the GAMA 
data, lTavloretd] ( l201lh obtained a = -0.68 and b = 0.73. This 
is close to the values ob tained from fitting to SDSS colours and the 
iKauffmann et alj j2003h stellar mass estimates, for example. 

Figure [T2l^a) shows GSMFs from GAMA data testing the ef- 
fect of varying colour-based M/L;. The values of a and b are such 
that the M/Li is 2.0 for galaxies with g ~ i = 1.2. The GSMFs 
show the flattening around and below lO^'^'^Mo, with a steepening 
below 10^-^ Mq: this is more pronounced with 6=1 than 6 = 0.5. 
Also shown is a comparison with the results of BGD08. There is 
generally good agreement between the GAMA and BGD08 results 
except at < lO^'^M©. This is despite the fact that the BGD08 re- 
sults are expected to be less complete in terms of SB. As noted 
above, this is because of the distance model used for the redshifts 
in BGD08. If the distances are changed to the model used here 
then there is good agreement. The GAMA results are more reli- 
able at ~ 1O*M0 because of the minimal dependence on the dis- 
tance model. For galaxies with 10* < M/Mq < 10*", 90% are 



Table 1. Galaxy stellar mass function. The values for masses lower than 
10*Mq should be regarded as lower limits (see ij |3.2) . The errors quoted 
are pseudo Poisson errors derived from the square root of the sum of weights 
squared. 
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brighter than Mr — —16, the GLF or GSMF is not significantly 
affected by uncertainties in the distances (Fig.[3j. 

Figure I12rb) shows the GSMF from the stellar masses of 
lTavloretal.l ( l201lh . strictly the fitted M/Li ratios in the AUTO aper- 
tures applied to the flux derived from the Sersic i-band fit (binned 
GSMF given in Table and the GSMF derived using the best- 
fit colour-based M/Li. These are nearly the same suggesting that 
a colour-based M/Li is easily sufficient for determining a GSMF 
assuming of course that it is calibrated correctly. From the GSMF, 
the total stellar mass density is 2.3 x lO^M© Mpc~^. This gives 
an fistars value of 0.0017 (relative to the critical density), or about 
4% of the baryon density, which is on the low-end of the range of 
estimates discussed by BGD08. 

3.4 Comparison with the A'-band galaxy luminosity function 

In order to compare with the shape of the GSMF, we also deter- 
mined the if-band GLF using the same V^ax values. For this, we 
used the A-band magnitude defined hy K = A"auto — iauto + 
zsersic, where the AUT O photometry is from the r-defined cata- 
logue ( lHiUetalJ|20Tll) . The reason for this definition is that for 
low-SB galaxies an aperture is more accurately defined in the SDSS 
r-band (or i-band) compared to the UKIDSS A'-band. This K — i 
colour is added to our fiducial i-band Sersic magnitude in order to 
be get a robust estimate of total A-band flux. The resulting GLF 
was simply converted to a GSMF using M/Lk = 0.5, which was 
chosen to give approximate agreement with the GSMF derived us- 
ing the iTavlor et alj ( I2OI ih stellar masses. The number densities 
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Figure 12. Galaxy stellar mass functions, (a): These GSMFs were derived from Petrosian magnitudes. The symbols with error bars represent the GAMA data 
using a color-based M/L^ (Eq.|5); the colours were derived from the GAMA AUTO photometry. The red and blue solid lines were computed using different a, b 
values. The red and blue dotted lines (offset xO.l) represent the results from BGD08 and the recomputation using the same sample after changing the flow 
model to that used here, respectively, (b): The symbols use the fitted GAMA stellar M/L^ from iTavlor et al.l J20T l|) applied to the i-band Sersic fluxes, while 
the dotted line uses a color-based M/L^. The solid line (off set xO.l) shows eff ectively the K-hanA GLF with K = JCauto — iauto + 'Scrsic and applying a 
constant M/Lx. while the crosses with error bars show the lPriver et al.l i2012l) ii'-band GLF. 



were divided by an average completeness o f 0.93 becau s e of th e 
reduced coverage in the A'-band [fig. 3 of ISaldrv et alJ ( l201(]h l. 
This scaled GLF is shown by the blue fine in Fig. llir b). Note that 
strictly the V^^ax values should be recomputed because of the dif- 
ferent coverage across the regions but this should have minimal 
impact on the shap e. We also show the GAMA iC'-band GLF from 
iDriver et alj ( I2OI2I) , which was derived from a different sample 
(0.013 < < 0.1, rpct < 19.4 and Kab < 18.1 with r-defined 
Ka,uto magnitudes) with the same M/Lk applied. 



The flattening from 10 "'^Mq to ~ 10 "Mq and upturn be- 
low these masses shown in the i-band derived GSMF is also seen 
directly in the iC'-b and GLF [Fig.[T2{b)]. Though in the case of the 
[Driver et alj ( 1201 2h result (standard Knax) it is less pronounced. 
This is an important confirmation of this upturn since, while there 
is some variation in M/Lk, the A'-band GLF is often used as a 
proxy for the GSMF. Previous measurements of the K-band field 
GLF had failed to detect this uptu rn using 2MASS photometry 
down to Lk < W^Lq jCole et al. 200 h Kochanek et al. 2001) 
or us ing UKIDSS with SPSS re dshifts l lSmith. Lovedav. & Crosi 
120091) : see fig 14. of ISmith et alj for a compilation. These measure- 
ments nominally probe far enough down the GSM F IO^Mq) 
that th e upturn should have been noted. We note thatl Merluzzi et al] 
( l2010h 's measurement of the if-band GLF in the z = 0.048 Shap- 
ley Cluster shows an upturn particularly in the lower-density envi- 
ronments, however, this does rely on statistical background subtrac- 
tion. The explanation for 2MASS-based GLFs missin g this could 
be the surface brightness limit. However, GAMA and lSmith et alj 
both used UKIDSS photometry. The difference in this case is that 
GAMA has redone the near-IR pho tometry using r-band defined 
matched apertures l lHilletal.ll20Tlh , and the magnitude limit is 
higher meaning the galaxies are typically further away (smaller on 
the sky) making near-IR photometry more reliable. 



3.5 The double Schechter function 

The shape of the GSMF is well fit with a double Schechter function 
with a singl e value for the break mass (M*), i.e. a five-parameter 
fit (BGD08: |Pozzetti et alj|2010l) . This is given by 



dM = e"^/' 




where dM is the number density of galaxies with mass be- 
tween A4 and M + dAi; with 02 < ai so that the second term 
dominates at the faintest magnitudes. Figure[T3]shows this function 
fitted to the GSMF data providing a good fit. The fit was obtained 
using a Levenberg-Marquardt algorithm on the binned GSMF be- 
tween 8.0 and 1 1.8 (Table[T]>, and the fit parameters are given in the 
plot. The fit to the lPozzetti etal. GSMF for z = 0.1-0.35 is also 
shown, which is similar. 

A natural explanation for this functional form was suggested 
by IPeng et al] (2010b). In their phenomenological model, star- 
forming (SF) galaxies have a near constant specific star-formation 
rate (SFR) that is a function of epoch. Then there are two prin- 
ciple processes that turn SF galaxies into red-sequence or pas- 
sive galaxies: 'mass quenching' and 'environmental quenching'. 
In the model, the probability of mass quenching is proportional 
to a galaxy's SFR (mass times the specific SFR). This naturally 
produces a (single) Schechter form for the GSMF of SF galaxies. 
Considering only mass quenching, the GSMF of passive galaxies 
is also determind to have a Schechter form with the same value of 
M* but with the faint-end (power-law) slope +1 compared to that 
of the SF galaxies. To see this consider a single Schechter func- 
tion GSMF and multiply by mass: M" M°'^^. Overall the 
GSMF of all galaxies is represented by a double Schechter func- 
tion with Qfi = Qf2 + 1- This is in agreement with our fit (Fig. 113b, 
which has ai — 02 ~ 1.12 ± 0.19. In fact, a good fit can be ob- 
tained by restricting qi = 122 + 1, making a four-parameter fit (at 
M > IO^Mq). 
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Figure 14. Color-magnitude diagrams. The points represent tlie GAMA 
sample used in this paper. The solid and dashed lines represent possible 
dividing lines between the red and b lue populations, wit h slopes of —0.03 
in (a) and using a tanh function from lBaldrv et al.l j2004) in (b). 



Figure 13. GSMF with a double Schechter fit to data at > 10" Mq. 
The data points represent the GAMA fitted stellar mass results. The solid 
line represents the fit to the data with extrapolation shown by the dashed 
line. The fit parameters are s hown with their l-sig ma errors. Also shown is 
a fit to zCOSMOS data from lPozzetti et alj j201(l . 

In the model, environmental quenching does not change the 
overall double Schechter shape as some SF galaxies are turned to 
red across all masses. The GSMF of the SF population remains 
nearly the same shape while the red-sequence GSMF has a scaled 
'copy' of the SF GSMF added so that it follows a double Schechter 
form most obviously in high density regions. 

To illustrate the origin of the double Schechter s hape of the 
GAMA GSMF as suggested by the lPeng et"al] ( l2010bl) model, we 
divided the galaxies into red and blue populations based on color- 
magnitude diagrams. Figure [T4l shows the g — r and u — r color- 
magnitude diagrams both versus Mr, with thre e possible divid- 
ing lines using a constant slope of —0.03 (e.g. Bell et al. 
and three using a tanh function (eq. 1 1 of Baldry et al. 2004), re- 
spectively. Figure [T5] shows the resulting red- and blue-population 
GSMFs with the dotted and dashed lines representing the six dif- 
ferent colour cuts (some extremely red objects were not included 
because the colour measurem ent was unrealistic , g — r > 1.0 or 
u — r > 3.2). Following the Peng et aEI ( 1201 Obi) model, we fit to 
the red and blue population GSMFs simultaneously with a double 
Schechter (ai, 02) and single Schechter function (a), respectively. 
The fits shown in Fig. [T5] are constrained to have: the same Ai* , 
02 = a, and qi = 02 -|- 1. A good fit with the five free parame- 
ters is obtained to the two populations when using a it — r divider. 
Note there is an excess of blue population galaxies above a sin- 
gle Schechter fit at high masses when using a g — r divider, the 
red population data were not fitted below 1O*'^M0 where there is 
significant uncertainty in the population type because of presum- 
ably large errors in the colours, and the inclusion of edge-on dusty 
disks is a problem for a sim ple red colour selection. Nevertheless 
the basic IPeng et alj ( 1201 O tf) model provides a remarkably simple 
explanation of the observed GSMF functional forms. 

3.6 The most numerous type of galaxies 

Are blue (irregulars, late-type spirals) or red (spheroidals, ellipti- 
cals) dwarf galaxies the most numerous type in the Universe (down 
to ~ IO^Mq)? Judging from Fig.[T5] the answer would appear to 
be the blue dwarf galaxies, i.e. star-forming galaxies. However, the 



i-i ^ (log 0*/lO"^ a) = (10,72. 0,71. -1,45) 




log („<t/Mo) + 2 log (Ho/70) 

Figure 15. GSMFs for the red and blue populations. The circles and squares 
with error bars, coloured according to the population, were derived using the 
divider that is shown as a solid line in Fig. lMf b). The solid lines represent 
fits to the data. The dotted and dashed fines represent the GSMFs using the 
six possible dividers based on g — r and u — r, respectively. 

measured number densities of both populations may be lower lim- 
its; and the measurement of the red population becomes less reli- 
able below about 10** '*Mq because of the smaller volume probed, 
the uncertainties in the colours, and the cosmic variance is larger 
because the galaxies are more clustered than the blue population. 

An alternative estimate of the number densities of red galax- 
ies can be obtained by considering the relative numbers of early- 
type galaxies in the Local Group, and then scaling the numbers to 
match the field GSMF at high masses (> lO^M©). This assumes 
that the Local Group represents an average environment in which 
these galaxies are located. Taking the catalogue of galaxies from 
Karachentsev et al. ( 2004), galaxies are selected within 1.4Mpc 
and with Galactic extinction less than 1.2 mag. The latter excludes 
two galaxies viewed near the Galactic plane (a biased direction in 
terms of detecting the lowest luminosity galaxies). The B-band lu- 
minosities are converted to stellar masses assuming: M/Lb = 3.0 
for early-type galaxies (RC3 type < 0); MJLb = 1.0 for M3 1, M33 
and the Milky Way, which have already been corrected for internal 
attenuation; and M/L_g = 0.5 for late-type galaxies (RC3 type > 6). 
From this, there are 6 galaxies with stellar mass > IO^Mq, which 
are M31, M32, M33, MHO, the Milky Way and LMC. For this 



© 201 1 RAS, MNRAS 000,[T]{T5] 



12 I.K. Baldry et al. 



0.100 



blue field galaxies • 
l.Ci early types (scaled) O 




8 9 
log (Jt/Mo) 



Figu re 16. GSMFs for the fi eld (GAMA) and the Local Group (derived 
from iKarachentsev et al.ll2004 ). The sohd Une shows the GAMA fit from 
Fig. [m with the dashed line representing the extrapolation. The dotted line 
represents the scaled LG GSMF in bins of 1 dex. The circles represent the 
blue field population while the diamonds represent the LG early types. 



population the Local Group, taken to cover a volume of lOMpc'^, 
is approximately 50 times higher density than the cosmic average. 

Figure [16] shows GSMFs for the field and the Local Group 
scaled to match, in particular comparing the blue field number den- 
sities with that inferred for the early types by scaling. It is likely 
that the LG sample is complete down to 1O^M0 wi th only some 
recently discovered satellites of M31, e.g. And XXI jMartin et alj 
|2009|) . suggesting that the bin shown here from 10^ to lO'^ is a 
lower limit. This analysis is consistent with the blue dwarf popu- 
lation being the most common galaxy down to lO^M©; at lower 
masses, it is not yet clear. 



3.7 Future work 

The GAMA GSMF is reliable down to lO^M© (corresponding to 
Mr ~ -16 with M/Lr = 0.5), which confirms the SDSS resuU 
(BGD08) with minor modification to the distances, assuming that 
the M/L values are approximately correct as a function of a galaxy's 
colour. In addition, there are ~ 350 galaxies in this GAMA sample 
between 10^ and IO^^Mq, and ~ 30 between 10^ and IO^Mq. 
There are a number of improvements to be made for the GAMA 
GSMF measurement at X < lO^M©: 

(i) The GAMA survey is ongoing with an aim to complete red- 
shifts to r < 19.8 over 300 deg^. This will approximately treble 
the volume surveyed for low-luminosity galaxies. 

(ii) There are about 2000 galaxies so far that have been spec- 
troscopically observed twice but with Q ^ 2. A careful coadd of 
the duplicate observations will yield additional redshifts for some 
of the low-SB galaxies. 

(iii) Flux measurements of currently identified low-mass galax- 
ies can be improved by careful selection of appropriate apertures. 
Automated Petrosian or Sersic fitting can lead to large errors for 
well-resolved irregular galaxies. 

(iv) Specialised searches can be made for low-SB galaxies that 
were missed by SDSS PHOTO, in particular, on deeper imaging 
provided by the Kilo-Degree Survey (KIDS) with the VLT Sur- 
vey Telescope and the VISTA Kilo-Degree Infrared Galaxy Survey 
(VIKING). In the longer term, a space-based half-sky survey, such 
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Figure 17. Size-mass relations. The blue and red solid lines show the 16th, 
median and 84th pe rcentiles of the effective radii from the GAMA Sersic fits 
teelvin et al.l201 l|) for the blue and red populations. The symbols represent 
measurements for irregulars and dwarf spheroidals. The dotted line outlines 
the region of incompleteness because of SB limits relating to the current 
GAMA sample. 



as that planned for the Euclid mission jLaureijs et allboioh of the 
European Space Agency, potentially will be able to detect low-SB 
galaxies with 10^ < M/Mq < lO" over 10^ Mpc^. 

The expected currently missed detection of low-SB galaxies 
is critical. In this sample, the observed number density for galax- 
ies with lO"'^ to lO'^Mo is only ~ 0.02 Mpc"^ dex"^ estimated 
using the den sity-co rrected l/^ax method. The predicted number 
bv lGuo et al] j201 ll) and by extrapolation of the double Schechter 
function is > 0.1 Mpc"'' dex~^. Thus we could be missing signif- 
icant numbers of larger low-mass galaxies. 

Figure [T7] shows the observed size-mass relation of galax- 
ies from GAMA for blue and red galaxy populations. For 
comparison, also shown ar e measurements of irregular galaxies 
jHunter & Elme gree i]|2006h using MfLy from t\\e B — V relation 
of Bell et al. ( 200J), an d Milky Way jOilmore et al.l2007h and M3 1 
dwarf spheroidals (e.g. iMc Connachie & Ir winl 20061 ; Martin et alj 
,2009.) using M/Lv = 2. The GAMA relation for the blue popu- 
lation follows an approximately linear relation above ~ lO^'^M© 
but appears to drop below the linear extrapolation at lower masses. 
The dotted line outlines the region where possible low-SB galax- 
ies missed by SDSS selection would be located. These would have 
Mr,50 ^ 24-25 mag arcsec^^ for the low M/L blue population. This 
is where an extrapolation of the mass-SB relation to low masses 
would lie (Fig.ll U. Thus it is essential to use a detection algorithm 
that is sensitive to these types of sources (e.g. Kniazev et al. 2005) 
at distances '--^10-50Mpc in order to test whether the lowest-mass 
bins are incomplete within the GAMA survey volume. For the star- 
forming population, obtaining redshifts is feasible but IFUs would 
be required if only part of each galaxy has detectable line emission. 

The low-redshift sample here only uses 5 per cent of the 
GAMA r-band limited main survey. The GAMA survey is also well 
placed to measure the evolution of the GSMF out to z ~ 0.6 for 
the most massive galaxies, study variations with environment and 
halo mass, and to study variations in properties with stellar mass. 
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4 SUMMARY AND CONCLUSIONS 

We present an investigation of tiie GSMF using tiie GAMA sur- 
vey. Throughout the paper, a recurring theme has been the ways 
in which different aspects of the analysis can affect the inferred 
shape and normahsation of a GLF or GSMF. In particular we have 
explored the importance of accounting for: bulk flows when esti- 
mating distances, large scale structure when estimating effective 
maximum volumes, the effect of using different photometric mea- 
sures, the surface brightness limit, and the effect of using different 
simple prescriptions to estimate stellar mass. 

The distance mo duli to apply to the magnitudes depend sig- 
nificantly on using the lTonrvetaUfeOOOl) flow model in compari- 
son with fixed frames (Figs. There is a noticeable effect on 
the measured number density of galaxies fainter than Mr = — 16 
(Fig. O. Using the same flow model with SDSS data brings into 
better agreement measurements of the GLF and GSMF between 
SDSS and GAMA [Fig.[9l;b), Fig.[T2l;a)]. For the same luminosity 
galaxies, the r < 19.4 (19.8 in G12) GAMA sample is less sensi- 
tive to whether the flow model is correct than the r < 17.8 SDSS 
sample. 

Measuring the GSMF accurately over a large mass range re- 
quires surveying a suitable volume > 10^ Mpc""^ to obtain at least 
tens of galaxies at the high-mass end (> lO^^M©), while the vol- 
ume over which low-mass galaxies are observed need not be so 
large. A problem arises in that the volume over which a galaxy is 
visible depends on its luminosity, and any variations in density as 
a function of redshift will distort the shape of a GSMF based on 
the standard Vinax method. Here we use a density-corrected K4ax 
method. This has been shown to be equivalent to a maximum- 
likelihood method I Colell201lh but is simpler to apply to an esti- 
mate of the GSMF. A volume-limited DDP sample of Mr < - 17.9 
(Fig. \^ was used to measure relative densities up to a given red- 
shift (Fig.|6}; and these are used to produce the density-corrected 
volumes (Eq. A useful diagnostic is to plot "K^ax versus Adr 
(Fig. |7}, which shows that K4ax increases nearly monotonically 
but with changes in slope compared to Vmax. The density-corrected 
V^ax method significantly reduces the difference in the measured 
GLFs between the regions compared to using the standard Viriax 
method (Fig. [Hi. 

There are small differences in the measured i-band GLF de- 
pending on the method of determining a galaxy's flux [Fig.|9ja)]- 
The AUTO apertures and Sersic fits recover more flux from early- 
type galaxies than the Petrosian aperture. This makes a significant 
difference at the bright end of the GLF. Converting the GLF to a 
GSMF using a colour-based M/Li relation results in a more obvi- 
ous flattening and rise from high to low masses as the b parameter is 
increased [Fig.[T2ja)]. Similar GSMF results are obtained whether 
using a fitted M/Lj for each galaxy or the colour-based M/L^ from 
lTavloretal.l l |201lh rFig. ll2r b)1. This is not surprising because the 
GSMF is only a one-dimensional distribution. We also find that the 
if -band produces a similar GSMF using a constant M/Lk = 0.5. 
This is an important verification of the upturn based on a simpler 
assumption that the K-hand approximately traces the stellar mass. 

As in BGD08 and IPozzetti et all feOlOh . we find that the 
double Schechter function provides a good fit to the data for 
A4 > W^Mq (Fig.ll3t. This is approximately the sum of a single 
Schechter function for the blue population and double Schechter 
function for the red population (Fig. [TS] ). This supports the em- 
pirical pict ure, quench i ng mo del, for the origin of the Schechter 
function bv lPeng et al] bOlObh . 

Blind redshift surveys, like GAMA, are better at characteris- 



ing the GSMF for the star-forming field population than the fainter 
and more clustered red population. In order to test whether the blue 
population is the most numerous in the mass range 10^ — lO^^M© as 
implied by the GAMA GSMF, we determined an approximate LG 
GSMF and scaled the resulting numbers to match the field GSMF 
at masses > lO^'M© (Fig. 1 161). The numbers of early types in the 
cosmic-average GSMF implied by this analysis are below that of 
the directly measured blue population. 

Accurately measuring the GSMF below lO^M© is key 
to probing new physics. For example, a simple prescription 
for preventing star formation in low-mass haloes, considering 
temperature-dependent accretion and supemovae feedback, results 
in a peak in the GSM F for star-forming galaxies at about lO^M© 
jMamon et al.ll20ll]) (note the overall baryonic mass function con- 
tinues to rise in their model). The problem with observing low-mass 
galaxies, 10^-10* Mq, is not the GAMA spectroscopic survey limit 
(r < 19.8) at least for the star-forming population but primarily the 
well-known issue with detecting low-SB galaxies (Fig.ll U. Thus a 
future aim for the GAMA survey is to characterise the extent of 
the missing ^ kpc size low-mass population (Fig. llVt . which ul- 
timately will require high quality deep imaging with specialised 
followup. 
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